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ABSTRACT 

Our concern is with the flow around an infinite cylinder, which at a certain instant is 
impulsively started to spin. The growth of vortices in the resulting boundary layer occurring 
outside the cylinder is investigated. This layer is essentially a Rayleigh layer which grows 
with time, so the mechanism involved is similar to that studied in Hall (1983). Vortices 
with wavenumber comparable to the layer thickness are shown to be described by partial 
differential equations. It is found necessary to solve the unsteady partial differential equations 
that govern the system numerically. We assume that the Rayleigh layer is thin so particles 
are confined to move in a path with radius of curvature the same as the cylinder. The Gortler 
number is a function of time, so we consider the time scale which produces an order one 
Gortler number. We consider the right hand branch calculation by letting the time tend to 
infinity, also inviscid Gortler modes are considered. 
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§1 Introduction 


The problem we are considering, is that of an infinite cylinder held fixed, until at a 
certain instant when a torsional impulsive velocity is imparted to the cylinder, without 
loss of generality this time may be considered to be t = 0. The cylinder is immersed 
in a fluid of kinematic viscosity v and density p, and the cylinder has radius a. The 
torsional velocity imparted to the cylinder is Q. It is known that a layer of thickness 
a(vt) 2 occurs at the cylinder, this is a Rayleigh layer which was shown to occur on 
a infinite fiat plate by Stokes (1855) and later by Rayleigh (1911). We consider small 
time scales so that the resulting layer is thin and thus the particles contained in this 
layer can be considered to move in paths with the same radius of curvature as the 
cylinder. We may define the Gortler number as 


G = 


Q 2 a 2 (vt) 2 
^2 


( 1 . 1 ) 


in the same way as in Seminara Sz Hall (1976) and Hall (1984). These papers considered 
centrifugal instabilities in Stokes layers, so (i 4) 2 is replaced by (u/u ) 2 , where u> is the 
frequency of the oscillations. The Reynolds number is defined here by 



( 1 . 2 ) 


The ensuing problem is similar initially to Hall (1983), who considered a spatially 
growing boundary layer occurring on a concave wall. We wish to investigate vortices 
similar to those found to occur by Taylor (1923) and Gortler (1940). In Hall (1982 
a,b) a self-consistent asymptotic description of the growth of short wavelength Gortler 
vortices was given. In particular in Hall (1982 a), it was shown that the equations gov- 
erning the centrifugal instabilities axe partial differential and not ordinary differential 
equations as discussed by Gortler (1940) and subsequent authors. In Hall (1982 b) the 
nonlinear development of the vortices was considered. In this paper we concentrate on 
disturbances with wavelengths initially of order one, but will consider the right hand 
branch in terms of local wavelength and local Gortler number. 


1 



We may now apply the usual non-dimensionalisation and non-dimensionalise time 
with respect to a 2 /i/, thus the Gortler number in (1.1) may be written in terms of the 
Reynolds number R e , and the non-dimensional time T as 

G = R\T~>, (1.3) 

To invoke the mechanism as described in Seminara & Hall (1976), it is necessary 

_ 4 

to regard G as order one, and this can be seen to occur if T ~ R e 3 , and since we 

consider R e 1, this is a small time scale. We may define a similarity variable r], and 

thus it is possible to solve the basic flow equations in the resulting layer, note that we 

_ 2 

consider a radial scaling of R e 3 so that 

where R is a dimensional radial variable. The essential point here is that the equations 
governing the motion in this layer involve temporal derivatives as well as radial ones, 
this is equivalent to ensuring ^-derivatives are retained in spatial problems, a fact 
that seems to have been first noted by Smith (1955). Smith (1955) retained some of 
the terms associated with the non-parallelism of the Gortler problem. For Tollmien- 
Schlichting waves Bouthier (1973) and Gaster (1974) incorporated extra terms to ap- 
preciate the effect of boundary layer growth, which greatly improved the comparison 
between experimental and theoretical results. Although it was noted by Gaster (1974) 
that the Bouthier- Gaster work just amounts to a successive approximation technique, 
Smith (1979) attempts to remedy this and does so effectively. Here we consider an 
unsteady instability, this is included in Kerczek & Davis (1974) and subsequently in 
Kerczek & Davis (1976) and Kerczek (1982), where in these papers they consider pe- 
riodic flow. The description of the flow in a boundary layer growing with time on a 
semi-infinite flat plate was attempted in Otto (1992 b), employing a triple deck struc- 
ture as in Smith (1979 a), so it is known theoretically that vortices occur in a time 
dependent growing boundary layer. 

Here we concentrate on disturbances, of order one wavelengths as mentioned pre- 
viously as these are most likely to be of practical importance. The development of 
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the vortices is caused by the inclusion of a small hump. This hump is considered to 
be the same size as the boundary layer, this technique was employed in Denier, Hall, 
& Seddougui (1991) in the spatially varying Gortler problem. Various forms of this 
hump are chosen and give rise to different regions of vortex activity, subsequently re- 
ferred to as vortex wedges. In this paper we do not have the freedom of Denier, Hall 
& Seddougui (1991) to introduce the disturbance at a certain spatial position as we 
have said that the impulsive motion starts at time equals zero, and the hump must be 
present for all time. 

The identification of the time station T, when a particular vortex wavenumber is 
neutral incurs the same problem as in non-parallel spatial calculations, refer to Hall 
(1983), Gaster (1974) and Smith (1979 a). In this paper the flow quantity chosen to 
be representative of the vortex size is 

OO 

E = Ju 2 dY, (1.5) 

o 

and thus the neutral position is defined by the position where a, defined by 


a = E - 1 


dE 
dT ’ 


( 1 . 6 ) 


becomes zero. 

Thus at this time station we may calculate the local wavenumber a and the local 
Gortler number G ? . The values of the equivalent pair in Denier, Hall, & Seddougui 
(1991) are found to depend on the initial form of the disturbance. The initial form of 
the disturbance is a function of the hump at the wall. 

We now go to investigate inviscid Gortler modes, where we have a? ~ 1 and 
Gf 1 described in Denier, Hall, & Seddougui (1991). We also consider the righthand 
branch problem such that for a.f > 1, we have Gf ~ a^,, as investigated in Hall (1982 
a). In this problem we have constant curvature, thus the right hand branch cannot be 
encountered by letting k ~ .V? as in Hall (1982 a). In order to investigate the right 
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hand branch it is found necessary to consider very small initial wavelengths A; at a 
time T such that, 

fc~T _5 with T 1. (1.7) 

We also consider matching between the inviscid solutions and right hand branch solu- 
tions as in Denier, Hall & Seddougui (1991). 

The approach adopted in the remainder of this paper is as follows; in section 2 we 
derive the linear 0(1) wavenumber regime problem, in section 3 the inviscid Gortler 
problem is tackled, section 4 contains details of the right hand branch calculations, 
section 5 includes description of the numerical techniques involved in the solution 
of the equations derived in sections 3 and 4. Finally section 6 contains results and 
discussions including comments on matching between the solutions of section 4 and 
section 5, and a discussion of possible further developments to take account of fully 
nonlinear right hand branch calculations. A discussion of extension to the problem to 
further geometries is also included. 
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§2 Formulation of linear order one problem 


We non-dimensionalise the velocity field with respect to ati, and the pressure 
terms with respect to pa 2 Cl 2 . The induced motion is taken to be independent of 0 , the 
coordinate system employed is that of cylindrical polars with components (r, 9 , z), and 
corresponding velocity components ( u,v,w ), the Navier-Stokes equations then become, 


J_dU dU_ w dU_ _V^_ _ _dP_ J_ 
R e dt + dr + dz r dr ^ R e 


(( d 2 d 2 \ id i \ __ . . 

VVSr 2 + dz* ) + r dr r 2 ) U ' 1 ^ 


1 dV TT dV dV VU 
T e ~dT + u ’fr +w 'K + ~ 




d 2 d 2 \ 


1 dw dw Tir dW dP 1 

■ 5 — 5T + ^ 1" ^ ~a — — ~~5 1" TT 

R e dt dr dz dz R e 


+ -§-~\) V ’ 

r dr r l J 

((d 2 d 2 \ 1 

\\dr 2 + dz 2 ) + r dr) W ' 


dr 2 + dz 2 J 
d 2 


dU_ U dW 
dr ^ r ^ dz 


(2.16) 

(2.1c) 

(2.1d) 


We now introduce an axi-symmetric hump at the cylinder, defined by r = 1 + 6f (z) 
. Following the arguments in the introduction, it is found necessary to look at a time 

_ 4 

scaling T ~ R e 3 so that the Gortler number is order one, and in this section it will be 
taken to be one with no loss of generality. Thus we define boundary layer variables, ij 
and Z given by 


r = 1 + R e 3 t], 


Z — Re 3 




(2.2a) 


and a time boundary layer variable r, defined by 

T = R7*t, 


(2.26) 


To determine the basic flow we assume it to be independent of Z and a function of £, 
a similarity variable given by 





(2.3) 


The basic flow is found to be governed by 


d 2 v n/ .dv 


= 0 , 


(2.4a) 
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with boundary conditions 

v = 0 as £ — > oo, 
v = 1 at £ = 0, 

The solution of (2.4) is the error function, given by 

OO 

5(0 = 7?/ e ~ f4 
£ 


We now introduce a disturbance so that the flow field is now given by 
u = (0, v, 0, p ) + 6(R e 3 u,v,R e 3 w,R e 3 p) H , 


(2.46) 

(2.4c) 


(2.5) 


( 2 . 6 ) 


where the quantities with tildes are functions of rj , Z and r, and, as previously stated 


independent of 6. We take 6 to be a small vanishing parameter, so that we are in a 
linear regime, with this assumption the governing equations are found to be, 


du dp 

2vv = + 

or or) 


d 2 d 2 \ 
dr) 2 + dZ 2 ) 


u , 


(2.7 a) 


dv dv „ / d 2 d 2 \ _ 

!h + !h l u - \d^' ir dz 2 ) v ' 


dw 

dr 


dZ ^ \dr) 2 dZ 2 ) 


du dw 

!hj + dz =0. 


w, 


with boundary conditions given by 


(2.76) 
(2-7c) 
(2.7 d) 


u = v = w = 0 as r) — > oo, 


so that the disturbance is confined to the Rayleigh layer, 

dv 

U = v = — = 0 at 1) = f(Z). 
Or) 


(2.7c) 


(2.7/) 


Here (2.7 f) represents the no-slip boundary conditions at the hump on the cylinder, 
which is chosen to be scaled with the boundary layer thickness multiplied by the dis- 
turbance size 6. Note that since we looking close to the cylinder the terms proportional 
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to secondary radial corrections are neglected, such as 1/r Jb. Let us now apply the 
Fourier transform to the disturbance quantities, thus the transform U of u is defined 

by 

OO 

= vb / ie ‘ kZdZ ■ < 2 - 8 > 


U 


Thus the system (2.7) now becomes 


= _i k p + (&__ k 2\ W , 


dr 


d r? 


^- + ikW = 0 , 


(2.9a) 

(2.96) 

(2.9c) 

(2.9d) 


with boundary conditions of 

U -V = W -i) 

U = W = 0, V = -F 

where F is the Fourier transform of the hump f(z). Note that is the absence of the 
hump there would be no vortex mechanism. Now by elementary manipulation, we 
may eliminate W and P, to obtain the coupled partial differential system, 


as r] —> oo, 

(2.9e) 

dv 

ft? at „ = 0. 

(2.9/) 


di f dr 


d 2 


a 2 

dr] 2 


-k 2 )U = 2k 2 vV, 


_ jfc2 _l\v = —U 

Or] 2 dr) dr] ’ 


(2.10a) 

(2.106) 


this system is the Gortler equations for a vortex with temporal growth, with Gortler 
number equal to one. It is now noted that there will be a discontinuity in the dis- 
turbance flow field due to the impulsive start, so it is necessary to solve the problem 
asymptotically for small r. Since the solution of (2.10) necessitates a numerical ap- 
proach, this asymptotic solution may be used as an initial condition. It can be seen 
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from (2.10 b) that V^G, r) ~ r 2, Z7(0,r ) ~ r 2 for t<1 following this scaling we 
seek a solution of the form 


v = T-?v ( b(0 + 


(2.11a) 


C/ = r2C7 0 (^) + ... , (2.116) 

where £ is the similarity variable defined by (2.3), it may be shown that Vq(£) satisfies 


a 2 Vo dV 0 

V + 2£ lf +2V " =0 ’ 


(2.12a) 


with boundary conditions 


Vb(0) = 2 


m 


and Vo ( 00 ) = 0, 


(2.126, c) 


where F(k ) is the Fourier transform of f(z), the form of the hump. The solution of 
(2.12) is given by 

Fo=2^e-<’, (2.13) 

v 71- 

Note that we take the radial component Uo to be zero for small r, as it is 0(t 2 ) smaller 
than Vo. The equations may now be solved numerically, the solution is discussed in 
section 5. Note that since this system is linear F(k ) may be taken to be unity, this 
corresponds to a 6-function shaped hump, so is relevant to a localised imperfection on 
the surface of the cylinder. 


§3 The inviscid limit of the Gortler problem 

In Denier, Hall & Seddougui (1991) the inviscid Gortler problem was discussed, 
it corresponds to a large Gortler number and order one wavenumbers. The temporal 
inviscid Gortler problem is virtually identical to the inviscid problem set as an exercise 
by Drazin & Reid (1979) and solved by Floryan (1986). Thus we have the governing 
equations (2.10) with a Gortler number other than one, taken to be G 

{w- k 2 -£)(w- kl ) 0=2dk2 ™’ (3 ' lo) 
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(3.16) 


(w- k 2 -l) V = T/’ 


We now consider a solution of (3.1) in the inviscid regime, that is G 1 with k = 0(1). 
It may be shown that the growth rate is proportional to following Denier, Hall & 
Seddougui (1991). Thus the disturbance in this regime takes the form 


(U,V) = E(U,G~?V), (3.2 a) 

where 

E = exp J ft H ^ , (3.26) 

This form of the disturbance may be substituted into (3.1), and the two equations can 
be combined to produce, 

?{w- k *) u - 2k ' u *%;= 0 ’ (33) 

In Denier, Hall & Seddougui (1991) the equivalent equation has a closed form of 

solution. Unfortunately (3.3) does not, so thus must be solved numerically. The 

equation (3.3) was solved using a fourth order Runge-Kutta scheme. The growth rate 

is plotted against k in figure 1. It is interesting to look at the solution of (3.3) as k 

increases to infinity, thus effectively encountering the right hand branch. We anticipate 

1 

that the growth rate tends to a constant /? 0 2 , so that we expand the growth rate as 

P = (Po + k * Pi H ^ (3-4), 

the size of the second order term is found from the size of the layer found to occur 
adajacent to rj = 0, from elementary observations of the numerical results. 

Thus we introduce a boundary layer variable x, given by 


1 


k-* x . 


(3.5) 


In this layer it is known that the basic flow v tends to 

1-A*“£x + -" > 


(3.6) 
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where A = — Vj,|^=o = Expand the disturbance U as 


U = U 0 + kslh + 


at leading order we obtain 


— /3 0 U 0 -f- 2,\Uo — 0. 


This is a compatibility condition, so for non-trivial solutions we require /3 0 = 2A, at 


next order we obtain 




For decaying solutions as x —■ y oo> the solution of (3.9) is given by 


Di-JMip (*+*■}. 


(3.10) 


To solve the boundary condition at x = 0 we require that 


(3.11) 


thus we can find the second order growth rate /?i in terms of the eigenvalues of the 
Airy function. Notice that this function has many zeros on the negative real axis, we 
are interested in the highest value of the growth rate, that is the first eigenvalue. Thus 
as k — + oo, we have 

P — > a/ 2A — ^pAH"3 4 , (3.12) 

3 3 


which we compare with the right hand branch growth rate in section 4, this result is 
also included in figure 1, note this agreement is reasonable, but it should be noted that 
k~a term does not tend to zero very quickly. Also included in figure 1 is the righthand 
growth rate value, and it can be seen that the solid curve is tending to this value. 


§4 Right hand branch calculation 

We consider the right hand branch where k is large, and large Gortler number, 
G as introduced in section 3. Neutral vortices in this regime were described by Hall 
(1982 a). In the neighbourhood of the neutral position, temporally amplified modes 
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with growth rate 0(k 2 ), are found to occur. In this regime, we know for neutral modes 
G ~ l’ 4 , so we write 

k = (4.1) 


We note that the vortex activity is likely to be concentrated where the Rayleigh cri- 
terion is most violated, for the problem considered here this occurs at the wall, which 
is the same as the Taylor problem, one of the problems considered by Hall (1982 a). 
We consider a disturbance proportional to 

E = exp J Po (<f>) + e Pi ( <t> ) H , (4.2) 

where 1 /e = k, and thus G may be expanded as 

OO 

G=Y,9n‘~ i+n - (4-3) 

n=0 


Thus the governing equations become 


d 2 


dr] 2 


dr 

dr} 2 


d 2 


^-1-^ U = 2e>GvV, 


dr} 2 

d 


dr 


dv 

v = u 7r . 

orj 


(4.4a) 


(4.46) 


Here we consider non-neutral solutions, thus /?o is non-zero it is possible to derive 
a compatibility relation, so in (4.4) we now have 


0 U,V) = (u 0 ,k- 2 v 0 )E , 


(4.5) 


We introduce a layer of thickness A* 2 G » , that is k 2 . So substituting (4.5) into 
(4.4) in this layer, located at the cylinder, we obtain 

(At 2 +/?)A t 2 u 0 = 2u 0 , (4.6a) 

(A t 2 +/?)u 0 =Au 0 A t2 (4.66) 

where we have assumed for rj <C 1, v — > 1 — A 77 + • • •. Thus combining (4.6 a) and (4.6 
b) we obtain, 

(A t2 + /?) - 2A = 0. (4.7) 
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Thus as A* — * 0, the growth rate j3 tends to \/2A, which agrees with the inviscid 
problem as k — » oo. Note that we also have an order one cut off point where the 
growth rate is zero given by 

A f = (2A)*. (4.8) 

Thus it is now possible to continue in the same manner as Hall (1982 a) and derive 
the form of the neutral modes, this analysis is not included as the conversion to this 
case is trivial and has the same conclusions. Note that in Denier, Hall & Seddougui 
(1991) they found that the position of the vortex activity varies with A*, we do not 
have that here because our growth rate is temporal and associated with d/dr, whereas 
in Denier, Hall & Seddougui (1991) the growth rate is associated with ud/dx. They 
also found as A* — ► 0, that there is a singularity in the growth rate and thus there is 
an intermediate G* wavenumber regime required to match with the inviscid problem 
but that does not occur here, as we have direct correlation between large- inviscid 
growth rate and small- right hand branch growth rate. 

§5 The numerical work 

Now let us discuss the numerical scheme employed to integrate the inhomogeneous 
system for U and V. The scheme is based on a finite difference scheme in the 77 -direction 
together with an Euler step in r. Thus for example (2.10 b) becomes 

v; + ( %^vl+vii 

At y h 2 

where i denotes grid position in the ? 7 -direction, whilst ± denotes the function evaluated 
at r + At and r respectively. Furthermore h and At are the step lengths in the rj 
and t directions respectively. If the disturbance velocity field is known at time t, then 
the tridiagonal system associated with (2.1) can be solved to give V at time t + At, 
a similar discretization of (2.10 a) then produces a pentadiagonal system for U , this 
system can be solved for U at time t 4- At; this procedure can then be repeated to 
march forward in time. This scheme is fully implicit so is stable for At = 0(h) rather 
than At = 0(h 2 ) which would be the case with an explicit scheme. 


Vih - V t r x 


- k 2 V~ + - 1 — 
* ' + At’ 
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After some grid checks we found that 300 steps in r/, with infinity r? = 20 was 
sufficient to present the given results to sufficient accuracy. The value of A r used 
to calculate the results given, was Ar = 0.001. The initial value of r where the 
calculations were begun was varied until it was sufficient small to calculate U and V 
with sufficient accuracy, the value used was Ar = 0.01. 

As an indication of vortex growth, the energy was computed, defined by 

OO 

E = J U 2 drj, (5.2) 

o 


as mentioned in the introduction and then the growth rate /? was calculated, defined 

by 

(5.3) 




For all our calculations /? was found initially to be negative so at sufficient value of 
time the flow is stable. At a later time, say tq , we found that /3 becomes positive until 
ri, where f3 again becomes negative. Thus there is a finite time interval where the 
basic unsteady boundary layer is unstable. At r = To and r = ri, the instantaneous 
wavenumber kr% can be plotted against the instantaneous Gortler number to 
generate the neutral curve for spin- up, this is shown in figure 2, for G = 1. 

We conclude that spin-up is instantaneously unstable for small time interval, 

_ 4 

~ R e 3 ; during this time turbulent flow will almost certainly be set up so that stable 
rigid body rotation may never be achieved in the practical situation. Furthermore 
we note that the induced velocity for any shaped bump can be found by numerical 
integration of the transformed velocity field. Several different humps were tried, and 
the results were found to be quantitatively similar, so we only include the plots for 
the humps included in Denier, Hall & Seddougui (1991). In the graphs the ordinate is 
the 2 -axis and the abscissa represents the time axis. Figure 3, represents the contours 
of maximum V component for the hump 7 r^ /4e~ z2 / 16 , and figure 4 the contours. for 
7 T 2 2 / 32 e -2:2 / 16 . Note these are basically the symmetric and anti-symmetric cases. 
The region of vortex activity is basically the same for any form of hump, even step 
functions. 
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§6 Conclusion 

We have shown for the spin-up problem, that Gortler vortices may occur when 
forced by an axi-symmetric hump at a cylinder started to move impulsively, also that 
the shape of the hump does not effect the amplitude of the vortices found to evolve, 
or greatly alter the time of the onset of the instability. 

We note that the growth rate, for an order one initial wavelength has two zeros, 
namely r = To and r = t\ , thus between To and t\ , /? the growth rate is positive and 
during this finite time interval spin up flow is unstable. The neutral curve is plotted 
for the humps, notice it is independent of hump shape as it is done before the Fourier 
inversion. We also produce figures showing vortex growth for each hump considered, 
the plots are of time against the maximum value of V for a given z, this quantity is 
taken to be representative of the vortex amplitude. Note that to graphical accuracy 
the wavelength of the vortices obtained corresponds to the critical wavelength from 
the relevant neutral curve. 

In section 3, we described the inviscid Gortler modes, which we found to have 
a growth rate of \/2A G*, as k — > oo, in section 4 the growth rate of the right hand 
branch in the limit k — » 0, was found to tend to 2A, thus we have matching of 
the growth rates, there is therefore no further structure brought on by a singularity in 
the growth rate, as was found to occur in Denier, Hall, & Seddougui (1991). 

In section 4 we did not include the description of neutral vortices but briefly we 
see that the first order disturbance term satisfies, 



()2 

- — V 0 - aV 0 
4 


= 0 , 


where 


and 



8 = t 4 x ( — 4 


Yi 


+ 


where 




V&JJ ’ 

V& + V&)) ’ 


( 6 . 1 ) 


(6.2a) 


(6.26, c) 
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where the superscript + denotes evaluation at r/ + , the position of vortex activity, in 
this case r} + equals zero. Thus Vo is given by a parabolic cylinder function, ensuring 
the required decay as x — * oo, if a = — | — m, thus we have 


go 



<7i=0 


9 2 


3c* 2 

2V 0 +V 0 + ’ 


(6.3) 


note the terms VJ^ are the coefficients of the Taylor series of the basic flow evaluated 
in a layer of thickness k ~? , situated at r] + . 

It is also possible to extend the small wavelength analysis to a fully non-linear 
calculation, performing a similar calculation to Hall & Lakin (1988). 

It should be also noted that we may extend this analysis to further geometries, 
the first obvious extension would be to the case involving a sphere, a brief outline of 
the method involved follows. Firstly it is necessary to invoke slightly different scaling 
for a sphere involved in spin-up about an axis of symmetry, the scaling are now 


w ~ W, u,v~aW 2 /u, 


~ a, t ~ u/a 2 


p ~ pW 2 


(6.4) 


These scaling are the same as those employed in the torsional case of Otto (1992 a), 
where to the frequency is now replaced with u/a 2 Thus the governing equations become 



du v du 
U dr r dd 



V 2 u — 


2 u 


d 


r 2 r 2 sin 0 dd 


(v sin 9) , (6.5a) 


vw _o / <-/ U U t/iy LLU 1 W* 1 Op I __9 2 dll 

— + R 2 A u — +~ — ) cot0 = — ^+1 v 2 u+-— - 

r 


dv 

di 


dv v dv uv 
d dr + r dd r 


1 dp 
rdd 


r 2 dO r 2 sin# 


, (6.56) 


dw 


dw v dw uw vw 


-^- + .R 2 U— + --^r + — + — Cot# = V 2 u;- 


w 


dt 


dr r d9 


r 2 sin 2 6 ) ’ 


]_d_ 2 1 d 

r 2 dr + r sin 6 dd 


(y sin#) = 0, 


(6.5c) 

(6.5d) 


where, 


V 2 


1 d ( 2 d\ 1 d ( . d_\ 

r 2 dr y dr J r 2 sin d dd y dd J 
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Note that we have assumed that the flow is independent of the azimuthal coordinate 
</>, as was done in the torsional case of Otto (1992 a), this also corresponds to assuming 
the flow in the cylinderical case is independent of 8. A Reynolds number has also been 
introduced, where this is given by Wajv. 

Prom similar arguments to those employed in the torsional case discussed in Otto 
(1992 a), the onset of a vortex stability is most likely to occur at the equator, thus we 
regard the basic velocity in this area. The derivation of the basic flow was attempted 
in Benton (1965), and he showed the basic flow components to be given by, 

w ~ sin 0/(r, t), (6.6a) 


v ~ cos 8 sin 8g(r, t). (6.66) 


Note in the region of 7r/2 we introduce a boundary layer variable 8 <C 1, and the basic 
flow quantities expand as 


w ~ f(r,t) + 0(8 2 )... , 

(6.7a) 

u ~ 0(8) + . . . 

(6.74) 

The Gortler number is given by 


COIN 

h 

cs <u 

II 

o 

(6.8) 


and so in a similar manner to the cylinder problem, we employ a time scaling so that 
the Gortler number remains of order one, thus we introduce 

r = l + R e 3 r?, t = R e a r, 8=— + R e 3 8. (6.9) 

Now in the region of the equator of the sphere the basic flow is found to be given 
by an error function, we perturb this basic state by the quantity 

3 ft, .Re 3 V, Re 3 , (6.10) 


As thus we arrive at the governing equations 

du _ _ dp ( d 2 d 2 \ _ 


(6.11a) 
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dv _ dp 
dr d9 


d 2 d 2 


dw _ dw 
dr U dp 


dr} 2 

d 2 


d9 2 

d 2 


dr} 2 + 89 2 


w. 


du dv 

— 1 r = 0 . 

dp d9 


(6.11 b) 
(6.11c) 
(6.11c?) 


Thus note if a Fourier transform is now applied to the these equations we arrive at the 
equations (2.9), with v replaced by w. Therefore we have shown that a sphere involved 
in spin-up is locally unstable to the same vortex structure as that found at a cylinder 
in the same situation. Although it should be noted that as in the torsional case of 
Otto (1992 a), difficulties may be experienced actually at the equator, so extension of 
the problem into the complex 0-plane may be necessary. 
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Figure 1: Inviscid growth rate curve, including asymptotic inviscid growth rate 
and the righthand branch growth rate value. 
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Figure 2 The 



Figure 2: Neutral curve of local GorfcL 
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Figure 4: Contours of the maximum of the V component of the vortex evolving 
with time, for an anti-symmetric hump. 
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